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Abstract 

Schramm-Loewner Evolutions (SLEs) have proved an efficient way to describe a single continuous 
random conformally invariant interface in a simply-connected planar domain; the admissible probability 
distributions are parameterized by a single positive parameter k. As shown in the coexistence of 
n interfaces in such a domain implies algebraic ("commutation") conditions. In the most interesting 
situations, the admissible laws on systems of n interfaces are parameterized by k and the solution of a 
particular (finite rank) holonomic system. 

The study of solutions of differential systems, in particular their global behaviour, often involves 
the use of integral representations. In the present article, we provide Euler integral representations 
for solutions of holonomic systems arising from SLE commutation. Applications to critical percolation 
(general crossing formulae), Loop-Erased Random Walks (direct derivation of Fomin's formulae in the 
scaling limit), and Uniform Spanning Trees are discussed. The connection with conformal restriction and 
Poissonized non-intersection for chordal SLEs is also studied. 



1 Introduction 

Critical systems in the plane (such as critical percolation, self-avoiding walk, ...) are generally conjectured to 
have scaling limits satisfying certain conformal invariance properties. Schramm-Loevifner Evolutions (SLE), 
introduced in ^23} . have proved a powerful tool to describe and study these limits. 

In several of these discrete models, such as percolation, one can define more than one interface. Scaling 
limits of these systems of interfaces should have similar conformal invariance properties. If one studies the 
limit of a single simple path, then conformal invariance and a "Markov" property specify the limit, if it 
exists, up to a positive parameter, denoted by k; this limit is SLE^. 

If one considers the joint law of several random curves, then the limit can be encoded through Loewner's 
equations as (one-dimensional) diffusions. The coefficients now depend on several parameters describing 
the boundary conditions. But these coefficients have to satisfy some compatibility conditions, expressed as 
commutation relations for associated infinitesimal generators. These commutation conditions can be recast 
as systems of linear PDEs with meromorphic coefficients and regular singularities along hyperplanes ([S!)- 
The key idea here is that we want to define distributions on systems of geometric paths, independently of 
any (artificial) time parameterization. 

A particularly interesting set-up is the following: (2n) points are marked on the boundary of a planar 
simply-connected domain, and these (2n) points are joined by n curves. Then it is shown in |S] that 
the possible systems of random curves can be obtained by solving a certain holonomic system, which is 
parameterized by n and n. 
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In this paper, we shall be mainly interested in the properties of this family of holonomic systems intro- 
duced in Solutions of this system can be related to particular (local) martingales of SLE. These local 
martingales are essentially the Girsanov densities for the systems of random curves w.r.t. independent SLEs. 
Our main goal here is to give integral representations of these solutions. When studying global properties of 
solutions of differential equations (asymptotics/boundary conditions, monodromy, . . . ), use of such integral 
representations is crucial. Also, although the theory of holonomic systems is rich and deep, getting solutions 
of a particular system is often challenging. 

Cardy's formula for critical percolation ([31123!) gives the probability that there exists a crossing from 
left to right in a rectangle. For k = 6, the system generalizes Cardy's formula to situations were sides of a 
(2n)-gon are set to alternate colors. Other related formulae have been proposed or proved, such as a formula 
for expected number of crossings (0), for crossing of annuli (0), Watts' formula for rectangles f^Tl IH]). 
and Pinson's formula for elliptic curves f|22j'). These formulae depend on one parameter (which is complex 
in Pinson's formula). 

For K G (0,8/3), the solutions of the system can be interpreted in the continuous limit in terms of non- 
intersection of n independent SLE^; and independent "loop-soups" f )16l I21[I29| ) with appropriate intensity. 
If K = 8/3, this loop-soup is empty, and we are studying the probability that n independent SLE8/3 do not 
intersect. As SLE8/3 is conjectured to be the scaling limit of Self-Avoiding Walks (^HDi one can think in 
terms of the scaling limit of non-intersecting self-avoiding walks. 

In the case k = 2, one can sec that the solutions of the system are exceptionally rational and connect 
them to Fomin's determinantal formulae for Loop-Erased Random Walks (see [201 El E])- This can also 
be seen as a particular case of the previous restriction construction (in the continuum). In the case k = 8, 
corresponding to the scaling limit of the Uniform Spanning Tree (UST, see |2Up. one can study the situation 
where alternate wired and free conditions on the boundary of a domain force the creation of n non-intersecting 
Peano paths on the Manhattan lattice. The scaling limit can be identified through Wilson's algorithm. 

It has been brought to our attention that the solutions obtained here appear in the so-called Coulomb 
Gas formalism (see [3, Chapter 9 in jH], and references therein), where they are referred to as Coulomb 
Gas representations. In this framework, deriving integral representations for correlation functions involves 
(chiral) vertex operators and screening operators, and the Feigin-Fuchs integral representation. Further 
developments along these lines pertain to monodromy, fusion, and BRST cohomology. 

In the present article, we obtain directly these integral representations starting from the associated 
holonomic system. The analysis of this system leads for example to results on rank and reductions using 
symmetry (as in Section 4.4). Also, we focus on the interplay between the algebraic, analytic and probabilistic 
aspects of the problem. For instance, delicate analytic questions regarding boundary values can be bypassed 
by probabilistic arguments (Section 6). In the case of Loop-Erased Random Walks and Uniform Spanning 
Trees (c ~ —2), direct connections with the combinatorial models are also established (3.2, 3.3). 

The set-up, objects and notations are essentially the same as in |S] and are mainly the standard SLE 
notations ((gt), k, . . . ), so we will recall them only briefly. This article can be read either as a self-contained 
study of a particular family of holonomic systems, or as a companion paper to 

The paper is organized as follows. After some background and notations, several situations leading to 
natural systems of n curves are described. Discrete models, such as critical percolation, loop-erased random 
walks (in relation with Fomin's formulae), and uniform spanning trees (in relation with Wilson's algorithm) 
provide important examples. In the continuous set-up, the theory of conformal restriction and loop-soups 
\'21\) can also be connected to the problem of commuting SLEs ([H])- In the following section, as a 
preparation, several particular cases (particular values of k and n) are briefly studied; in particular, for 
K = 6 and 2n = 6, we get new crossing probability expressions for critical percolation in conformal hexagons. 
Building on these particular solutions, general (formal) solutions arc given, as Euler integral representations. 
We then discuss the problem of identifying solutions corresponding to specific geometric configurations; or, 
in other terms, how to choose the cycle of integration in the Euler integral. 
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2 Background and notations 



First wc recall some definitions and fix notations. We shall be mainly interested in chordal versions of SLE, 
i.e. random curves connecting two boundary points of a plane simply-connected domain, satisfying some 
conformal invariance properties. For general background on SLE, see \2'6\ I28LITS] . Also, we will use at some 
point results on the restriction property and the "loop-soup" (see [161 1211 0^ ). 

Consider the family of ODEs, indexed by z in the upper half-plane H = {z : > 0}: 

o 

dtgtiz) 



9t{z)-Wt 

with initial conditions g{i{z) ~ z, where Wt is some real- valued (continuous) function. These chordal Loewner 
equations are defined up to explosion time (maybe infinite). Define: Kt = {z G M : < t}. Then {Kt)t>o 
is an increasing family of compact subsets of H; moreover, gt is the unique conformal equivalence ]H[\il't H 
such that (hydrodynamic normalization at oo): 

gt{z) = z + o(l). 

For any compact subset of H such that H \ iC is simply-connected, we denote by (j> the unique conformal 
equivalence H — > EI \ /iT with hydrodynamic normalization at oo; so that gt = 4>Kf The coefficient of 1/z 
in the Laurent expansion of gt at oo is by definition the half-plane capacity of Kt at infinity; this capacity 
equals (2t). 

If Wt = X -\- ^/nBt where (Bt) is a standard Brownian Motion, then the Loewner chain {Kt) (or the 
family (gt)) defines the chordal Schramm-Loewner Evolution with parameter k in (H, a:, oo). The chain Kt 
is generated by the trace 7, a continuous process taking values in H, in the following sense: H \ il't is the 
unbounded connected component of H \ 7[o,t] ■ 

The trace is a continuous non self-traversing curve. It is a.s. simple if k < 4 and a.s. space-filling if 
K > 8. Its Hausdorff dimension is a.s. 1 + k/8 if k < 8 (and 2 otherwise). 

If {D,x,y) is a simply-connected domain with two marked points on the boundary, SLE„ from x to y 
in D is defined as the image of SLE^ from to 00 in H (as defined above) by a conformal equivalence 
(H, 0, 00) {D,x,y) (that exists by Riemann's mapping theorem). With this definition, SLE satisfies a 
"Domain Markov" property, which, together with conformal equivalence, essentially characterizes it. 

In [5j , the question of defining several SLE strands simultaneously in a domain is addressed (say each of 
those strands is absolutely continuous w.r.t. SLE^). The key point is that an appropriate "Domain Markov" 
condition for the joint law imposes stark "commutation conditions" on the drift terms of the driving processes. 
Elucidating those conditions is the main object of We summarize the result in the following case: (2n) 
points xi, . . . , X2n are marked on the boundary of H, and we want to define jointly n SLEs connecting these 
(2n) points, such that the joint distribution depends only on the position of the marked points, is Mobius 
invariant and satisfies the appropriate Markov property. Then necessarily, the SLE growing at Xi is driven 
by: 

dxf ) = V^dsf^ + «;^(5,(xi), . . . . . .,gt{x2n))dt 

where i?'*-' is a standard Brownian motion and tp is a, non-vanishing function annihilated by the operators: 

J2k - n{l - 6/k) 

J2k ^l9k - {1-6/k){xi-\ h X2n) 
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From a CFT point of view, this corresponds to differential equations for correlation functions of (f>i.2 
primary fields. 

If n = 1, ip{xi,X2) = {x2 — xi)^^^/'', and: 

dxf ^ = V^dB[^'> + -^j^^ dt 

which is ordinary chordal SLE^ from xi to X2 (with a homographic change of coordinate). 

The main question addressed in this article is how to obtain general integral representations for solutions 
of this system and their interpretation in terms of scaling limits for different discrete models. 

3 Probabilistic interpretations 

In this section, we discuss probabilistic situations giving rise to natural examples of systems of geometric 
paths. We start with critical percolation, and then consider two other discrete models with particular 
harmonic properties, the Loop-Erased Random Walk (LERW) and the Uniform Spanning Tree (UST). We 
then proceed to show how to use the Poissonian structure of the loop-soup to define a natural notion of 
non-intersection for n chordal SLE^^'s, where < k < 8/3. 

3.1 Crossing events for critical percolation 

We discuss here consequences of the locality property for SLEg / critical percolation, in particular regarding 
holonomic systems. 

First, we describe a family of percolation events that appear as particularly well-suited for an SLE 
analysis. Recall the set-up of Cardy's formula, relating to critical percolation in a conformal quadrilateral. 
For simplicity, we will consider site percolation on the triangular lattice: each site is colored in blue or yellow 
with probability 1/2, all sites being independent. (Alternatively, the hexagons of a honeycomb tiling are 
colored in blue, yellow). A subgraph of the triangular lattice, with mesh e \ 0, approximates a quadrilateral. 
As boundary conditions, a pair of opposite edges of the quadrilateral is set to blue, while the other pair is 
set to yellow. For (plane) topological reasons, either the two yellow edges are connected by a yellow path on 
the lattice, or the two blue edges are connected by a blue path (see Figure 0. 

It seems quite natural to generalize this construction to (2n)-gons, n > 2. So consider a bounded Jordan 
domain D C C; ai . . . a2n are distinct points on the boundary (in counterclockwise order, say). The boundary 
conditions are alternate: for 1 < i < rt, the arc (a2i,a2i+i) is set to yellow, and (02^-1,02^) is set to blue 
(with cyclical indexing, i.e. a2n+i =0-1)- We are interested in the connectivity properties of the random 
percolation graph that are observable from the boundary. Denote by the edge (a,;, a^+i). Let c(ei, ej) = 1 
if ei,ej are of the same color and are connected by a path of this color, and c{ei,ej) = otherwise. Note 
that a path can include a portion of the boundary. An elementary event is an event of type 

C(e) = n {c(ez,ej) = e,,j} 

l<i<j<2n 

where e = (fi j")i<i<j<2n, ^i.j S {0, 1}. We now describe the non-empty elementary events. We suppose that 
the mesh e is small enough, so that the edges are nonempty and disjoint. From the Russo-Scymour- Welsh 
theory (and FKG inequality, see [TJ ) , the elementary events that are topologically possible will happen with 
probability bounded away from (and from 1) as the mesh e goes to zero. 

Let n > 2. There are exactly C„ (non-empty) elementary events with probability bounded away from 
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Figure 1: Two interfaces (red, fuchsia) bound a crossing of the lozenge 



as e \ 0, where C„ denotes the n-th Catalan's number: 

V n / 

n + 1 

Indeed, it is easy to see that there is a one-to-one correspondence between elementary events and non-crossing 
partitions of the set of blue edges. As is well known, this number is C„. This can be seen as follows: consider 
the smallest i > 1 such that (a2i-i,a2i) is connected to (01,02) by a blue path; this leads directly to the 
recurrence relation for Catalan's numbers. For instance, for n = 3, one gets C3 = 5 possible configurations 
for a conformal hexagon with alternate boundary conditions, as illustrated by Figure [3 

Alternatively, at each point where boundary conditions change, one can start an exploration path 
f |25[ 1261 ) that winds between the connected components of (ai_i,ai) and (0^,0^+1). Such a path is simple 
and ends at Oj, for some j [i and j have opposite parity). So the {2n) boundary points are paired by n 
non-intersecting exploration processes. This pairing is random; the number of such pairings (satisfying the 
non-crossing condition) is C„. Note also that a site that touches two exploration processes is pivotal for 
these events. For bond percolation on Z^, it is convenient to represent the interfaces between connected 
components in and connected components for the associated percolation configuration on the dual graph. 
The collection of these interfaces (with appropriate boundary conditions) is made of closed loops and n 
non-intersecting simple paths connecting the (2n) "free" points on the boundary (see Figure 

The exploration process for critical site percolation on the triangular lattice converges to SLEg in the 
scaling limit (HHUD)- It is conjectured, and supported by numerical evidence, that it is also the case for 
critical bond percolation on (and also for more general lattices). 

From conformal invariance of the scaling limit of critical percolation, the probability of any of these 
elementary events should define a function on the corresponding moduli space (i.e. the space of Jordan 
domains with (2n) (distinct) marked boundary points modulo conformal equivalence). Denote by M2n this 
moduli space, which can be seen as a smooth (2n — 3)-dimensional manifold, for n > 2 (for a discussion 
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Figure 2: The five configurations for a conformal hexagon (schematic) 
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of SLE and moduli space, see e.g. JSl)- Considering all the possible configurations, one defines a function 
A^2n M*^". Assuming that this function is smooth, the dimension of the image is at most (2n — 3), which 
implies the existence of a large number of smooth relations between probabilities of elementary events. 
Indeed, the dimension of the moduli space. (2n — 3), is negligible compared with C„: 

(J _ Cn) 1 A/2^(2rt)^"e-^" /Y 

n + 1 n' (V27rn(ri")e-")2 VTr'n^/^ 

as follows from Stirling's formula. The sum of all probabilities is 1; the nature of other smooth relations 
between these probabilities is unclear. Letting a given edge of the conformal (27i)-gon shrink to 0, one 
gets a conformal {2n — 2)-gon, so that the corresponding part of the boundary of can be identified 

with A^2n-2- Considering this operation for any edge, it appears that any affine relation between the C„ 
probabilities is proportional to the trivial (normalization) relation. 

Assuming conformal invariance, one can give a differential characterization of these probabilities. More 
specifically, let f(xi,...X2n) be the probability (in the scaling limit) of any one of the C„ elementary 
events associated with the configuration (H,xi . . ■X2n), xi < ■ ■ ■ < x^n- Let i G {1 . . .2n}; we consider 
an infinitesimal percolation hull at xi. The boundary changes color at Xi\ Smirnov's key result (see 1221) 
is that the percolation exploration process started from xi converges to SLEg started at Xi (be it chordal 
SLEg to another boundary point, or radial SLEg to an inner point, since all these are equivalent for short 
enough times). Let (7ti)o<it<77 be the exploration process, and (ii'u)o<u<r; be the associated family of hulls, 
where 77 is such that the vertices Xj, j ^ i, are not disconnected at time 77. Then the elementary event holds 
for (H, xi, . . . , X2n) if and only if it holds for the (random) conformal (2n)-gon (H \ K^, Xi, . . . , "fu, ■ ■ ■ X2n), 
for any u € [0, 77]; assuming conformal invariance, this has probability f{gu{xi), . . . , Wu, ■ ■ ■ guix2n)), where 
((7„) are the conformal equivalences defining the SLEg process. So {f{gu{xi), ■ ■ ■ ,Wu, ■ ■ . 5u(2^2n)))o<«<)) is 
a martingale. To sum up, assuming that the probability of an elementary event defines a smooth function 
on A^2)i, then the function / is annihilated by the following differential ideal: 

/„ = (/:i,...£2„,^-i,4,^i) 

where Ci = idu + x -x- ^J" (infinitesimal generator for SLEg growing at Xi), and £k = ^X^^i^'^^i- 

Note that the (real) Lie algebra generated by the {ik)k£Z is (isomorphic to the) Witt algebra, and that 
the subalgebra {£-i,£o,ii) is isomorphic to s[2(M), the tangent algebra of the Moebius group (the group 
of conformal automorphisms of a simply connected domain). The fact that one can explore the (2n)-gon 
starting at any of its vertices is a feature of locality. 

It is readily seen that /„, which is an ideal of differential operators with rational coefficients, is holonomic, 
i.e. the vector space of functions annihilated by /„ in the neighbourhood of a generic point has finite 
dimension (the rank of /„ is the dimension of this space). For background on holonomic systems, see e.g. 
|3H . It is elementary that the rank of /„ is no greater that 4". Indeed, if /, defined in a neighbourhood U 
of a generic point, is annihilated by /„, let F : [/ — > R^"'^^ , where: 

F=(9r ...9^^"/).e{o,i}-- 

Using the operators Ci, it readily follows that one can write diF = MiF for i = 1 . . . 2n, where Mi is 
a matrix of rational coefficients. A local solution of this system in a neighbourhood of {xi, . . . ,X2n) is 
entirely determined by the initial condition F{xi, . . . X2n)^ so the solution space is of dimension at most 
2^" = 4" (it would be exactly of dimension 4" if the Frobenius integrability conditions were satisfied: 
djMi - diMj = MjMi - M^Mj). Since F(o,...,o) = /, the rank of /„ is at most 4". Note that we have not 
used the operators £1, so this is a crude estimate. Alternatively, it is easily seen that the dimension 

of the characteristic variety of /„ is (2n), implying holonomy (see e.g. |U). Making use of conformal 
invariance, if / is annihilated by /„, then one can write: 

/(Xi, . . .,X2n) = g{[xi,X2,X3,X4], . . . , [xi , X2 , Is , a;2„] ) 
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where g is a function of the (2n — 3) variables j/,; = [xi, X2,X'i^ where [.,.,.,.] is the cross-ratio. Since 

Ci, . . . C2n annihilate /, g is annihilated by operators Ci, . . . Cm-'i such that the only second order differential 
term featuring in Li is /dyj. Reasoning as above, this implies that the rank of /„ is at most 2^""'^. We 
expect the rank to be exactly C„. 

In the case n = 2, if a function / is annihilated by J2, then by conformal invariancc {i-if ~ £o,f = ^if = 
0), one can write 

f{xi,...,X4) ^ g{[xi,...,X4]) 
where [.,.,.,.] is the cross-ratio. Then / annihilates I2 if and only if g satisfies the second order ODE: 

3u(l - u)g"{u) + (2 - 4:u)g'{u) = 

which is a Fuchsian differential equation with three singular regular points on the Ricmann sphere (namely 
0, 1, 00), so it is essentially equivalent to a hypergeometric differential equation (see j^J, p. 26). It is easily 
seen that the rank of I2 is exactly 2; at a generic point, the solution space is generated by a constant 
function and by the function appearing in Cardy's formula. It is worth remarking that we could have chosen 
any cross-ratio (there are six of them) to get the same solution space, which translates into exceptional 
invariance properties of the considered ODE (under homographies that fixate the singular locus {0, l,oo}). 
In fact, this is the only second-order Fuchsian ODE (with singular locus {0, l,oo}) that is invariant under 
the substitutions u i— > (1 — u), w 1— > so that the three singular points play exactly symmetric roles. 

3.2 Loop-Erased Random Walks and Fomin's formulae 

In J Fomin considers a loop-erased version of the well-known Karlin-McGregor Formula. Problems per- 
taining to the scaling limits of Fomin's formulae are studied in ^B]- From Wilson's algorithm (|2ni)j we know 
that there is an exact identity between branches of a uniform spanning tree (UST) and loop-erased random 
walks (LERW). It turns out that the natural way to define non-intersecting LERWs is to consider disjoint 
branches of an ambient UST. 

More precisely, consider the following situation. A simply connected domain D of the plane, with, say, 
smooth Jordan boundary, is approximated by a subgraph of a lattice with mesh e. Then take a uniform 
spanning tree of this graph with wired boundary conditions. Consider now xi, . . . Xn n points on the boundary 
and j/i, . . . , ?/„ n points one lattice spacing away from the boundary, such that xi, . . . ,Xn,yn, ■ ■ - yi arc in 
counterclockwise order. We condition on the event that the minimal subtree containing ?/i,...,y„ and 
the boundary has no triple point in the bulk, and that the branch of yi connects to the boundary at Xi, 
i = 1, . . . ,n. Then take the scaling limit of these discrete paths. We get n non-intersecting paths connecting 
Xi to yi, i = 1, ... , n; each of them has density w.r.t. chordal SLE2 ( |2U) ) . 

In the discrete setting, let H be the harmonic measure for simple random walk on the lattice killed when 
it reaches the boundary. Then the probability of the event considered above is given by Fomin's formula: 

det{H{y,,{xj})) 

As follows from in the scaling limit (say in the upper half-plane H), the Girsanov density of the system 
of paths w.r.t. independent chordal SLE2(a;i yi) is given by the scaling limit of this determinant divided 
by the product of its diagonal terms: 




In the continuum limit, we will see later an extension of this situation to values of n between and 8/3. 
Note that this discrete construction is reversible (for slightly different boundary conditions). This situation 
further illustrates the close relationship between LERW and discrete harmonic measure and its application to 
scaling limits. In the next section, we shall make use of the connection between UST and reflected harmonic 
measure. 
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3.3 Uniform Spanning Trees 



If wc specialize the results wc shall obtain later in this paper to k = 8, the formulae become symmetric 
in the (2n) variables and have a nice Riemann surface interpretation. On the other hand, chordal SLEg is 
known to be the scaling limit of (the Peano curve of) the Uniform Spanning Tree (see [201 )• It is not hard 
to think of boundary conditions that enable to define multiple exploration processes in the discrete setting. 
From Wilson's algorithm ((SODi it is then possible to extract the drift terms. Though, the connection with 
our formulae is not quite immediate, so we shall give some details here. 

Consider a discrete lattice approximation of a bounded simply connected with, say, piecewise smooth 
boundary and (2n) marked points xi, . . . X2n on the boundary. The n boundary arcs (x2, X3), {xa, xs), . . . , (x2„, xi) 
are wired, and connected by an external wiring; the n other boundary arcs are free (i.e. wired for the dual 
graph). One then samples uniformly from spanning trees with these conditions. If we erase the external 
wiring, one gets n trees (corresponding to the n wired boundary components). Each vertex z is connected 
to one of the wired boundary component with probability given by the (graph) harmonic measure, from 
Wilson's algorithm. In the scaling limit, this should converge to harmonic measure with normal reflection 
on the free parts (see Figure 0}. 




Figure 4: Tree (solid), dual tree (dashed), Peano paths (curved) 

With these conditions, the Peano path starting at X2i ends at X2i+i , so that the boundary conditions fix 
a pairing of the (2?i) points; note also that this pairing is different from the one most natural for Fomin's 
formulae. Assume now for simplicity that the domain in the upper half-plane H, and xi < a;2 < ■ • • < X2ri- 
We want to prove that this situation corresponds in the scaling limit to SLEg's whose drift terms are given 
by log-derivatives of the "partition function" : 

l<i<j<2n -"^ l<i<j<n ''=^ [UjlliUi - Xj)] 

where C is a Cartesian product of Ci, . . . , C„_i, where Ci is a cycle circling clockwise round the segment 
(x2i-\iX2i). Note that one can choose a determination of the integrand on such cycles. In this expression. 
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rewrite the product ni<j(% ~ "0 ^ Vandermonde determinant. Define also: 

u^^^du I (u^xiY^^du 

— uj- = 

1/2' « / V 1/2 




for i € {1, . . .n — 1}. It is well-known that (wi, . . . , ujn-i) is a basis of abelian differentials of the first kind 
for the hyperelliptic curve of genus g = n — 1: 

= [S - Xi) ■ ■ ■ {S - X2n) 

and (Ci, . . . , Cn-i) is half of a canonical homology basis for this curve (see e.g. ^Jj for background on 
compact Riemann surfaces). Then can be written as a determinant of periods: 

V'(x)= H {x, - x.y/^et I f u;A= [] {x^ ~ x,)^^Uct 

l<i<j<2n j l<i<j<2n 

Consider the Peano exploration starting from xi and ending at X2n- We want to recover the form of ip 
from discrete arguments. First, we note that the situation can be seen as a degenerate case of a UST in a 
multiply-connected domain. Consider n — 1 compact holes Ki, . . . Kn-i in H, with the following conditions: 
two points Xi and X2n are marked on M; {xiX2n) is free and {x2nXi) is wired; the boundaries of the holes are 
wired, and all the wired parts are considered as wired together. One recovers the previous situation is the 
holes are segments close to R. If z is any bulk point, then Wilson's algorithm provides n — 1 martingales for 
the Peano exploration process. More precisely, let Ht be the remaining domain after time t of the exploration 
(for some time parameterization); and Harm is the harmonic measure in Hf with normal reflection on the 
free part of the boundary {{xiX2n))- Then: 

t ^ Mt — *(IIarm(z, dKi), . . . Harm(z, dKn-i)) 

is a vector- valued martingale (see [201 )• Following [20], and the work of Makarov and Zhan for SLE in 
multiply connected domains (see [12 )j one determines from this the driving process of the scaling limit. We 
use here the notations and conventions of |H]. So {gt) is a family of conformal equivalences Ht —^ M extending 
through the holes and with hydrodynamic normalization at infinity: 

dtgt = '^——,dXt = ^/ndBt + n-^dt 

gt - Xt V 

for some covariant function t/i of the configuration, which we want to identify. Let z be a point on the free 
part of the boundary close to x ~ xi] expanding at z = a;, one gets: 

Mt = (z - xfl'^at + {z~ xf'^bt + ■■■ 

the half-integer exponents coming from the normal reflection. The martingale condition at t = gives: 

+ + + ({z-xf/'-a+(z~xf'%+---)=0 

where a and b are now (vector-valued) functions of the initial configuration. Considering the terms in 
(z — x)^'^/'^ and (z — x)^^/^ in this equation, one gets the necessary conditions k = 8 and: 

+ Ad^a = 6b. 

So the drift term dx log V' can be recovered from a, b, that are coefficients in the expansion of the harmonic 
measure. Note that dx^ogip is a real number, and the condition is a relation on vectors. Now we want to 
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prove that -0 as defined above satisfies this equation in the degenerate case where Ki = (X2X3), . . . Kn-i = 
{x2n-2X2n-i) ■ We assume that n > 2 (the case n = 2 being comparatively trivial). 

Let P be the matrix: 



P = 



l<i,j<n 



and Q{z) be the vector ^{t^'i, ■ ■ ■ , Then it is easy to see that '^{P~^Q) is the image of the harmonic 

measure vector Hq by a fixed triangular matrix. In fact, the imaginary part will also give a martingale. So 
we can replace the harmonic measure vector by P~^Q in what follows, for simplicity. Hence, if (ei, . . . e„_i) 



is the standard basis of ' 



one gets: 



/2 



a = p-\2e,) l[{x. - x)-'/^b = p-' ( ( ^ ) e, + ^e^] - x)-' 

i>l \ \i>l / i>l 

where P~^Q{z) = {z — x)^/'^a + {z — x^^'^b + ■ ■ ■ . Using the fact that: 

det(P) = det(P) TT:{p-^d^P) 

and Tp ~ Y\i<i<j{^j ~ ^iY^'^ det(P), after simplifications, we have just to check that: 

8Tr(F-i9^P)ei +8P9:^p-^ei = 4e2. 

This follows from Pd^{P-^) + {d^P)P-^ = and the fact that d^Lo[^^ = (1/2 - i)Lo[ if i e {1, . . . , n - 2}, 
so that: 

( * * * \ 

-i 



d^P 







V 







P 



/ 



and the previous identity can be read from the first column of {dxP)P ^■ 



3.4 Non-intersection and the restriction property 

As pointed out in [Q, the restriction property for SLEg/3 can be put to use to get natural constructions of 
non-intersecting SLEs. For simplicity, consider first a simply-connected domain in C with four points marked 
on the boundary, say (H, xi, X2,X3, X4), where H is the upper half-plane. Consider two independent SLEg/3's 
in H, from xi to X2 and X3 to X4 rcsp., with traces 7,7'. Say (gt) is the family of conformal equivalences 
associated with the first one (time parameterization is unimportant here). The restriction property for the 
second SLE implies that the law of 54(7') conditionally on {7 n 7' = 0} is the law of an SLEg/s from 
gt(xz) to gt{xA) conditioned not to intersect gt{l(t,oo))- Using at the same time the Markov property and the 
restriction property for the first SLE, one gets that {gt{l),gt{l')) conditionally on {7(^7' = 0} is distributed 
as (7,7'); where 7 and 7' are independent SLEg/3's going from gt{jt) to gt{x2) and from gtix^) to gtixA) 
respectively, conditioned not to intersect. Indeed: 

>C(5t(7(t,oo)),fft(7')l7n7' = 0) = /:(gt(7(t,oo)),5t(7')l7(o,t) ny = 0, 5t(7(t,oo)) n 34(7') = 0) 

= /:(7,7'l7n7' = 0) 

Moreover, the system (7,7') has also a restriction property. This can be generalized to n SLEs connecting 
(2n) points on the boundary, and conditioned on no intersection between any two of them. This induces 
a topological constraint on the order of the marked points (aii, . . . , X2n) on the boundary. If the points 
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(cci, . . . , X2n) are in cyclical order on the boundary, one is interested in "non-crossing pairings" of these 
points (deciding that an SLE connects .t,; and Xj pairs Xi and Xj). The number of such pairings is given by 
Catalan's number C„. The value of the probability of non-intersection associated with such a pairing can 
be seen as a function on the moduli space A^2n, and satisfies (2n) evolution equations ([H])- 

If K € (0, 8/3), this can be extended using loop-soups of intensity = (6 — k)(8 — 3k)/2k (see [TCI 1^ 1. 
More precisely, consider (2n) marked points on the boundary, n independent SLE^ connecting these points, 
and independent loop-soups L2, . . . , Ln with intensity A^. If the pairing is non-crossing, there is a positive 
probability that no loop in Lj intersects Ui<j7' and 7-', j = 2, . . . , n. Conditioning on this event (and 
assuming reversibility), one gets a Markov property as above for each of the boundary points. Here one uses 
the Markov property of SLE, the restriction property for the loop-soup, the identity connecting loop-soup, 
SLE, and restriction measures (EH]), and the Poissonian nature of the loop-soup. Note that one can rephrase 
the conditioning in a symmetric fashion. Indeed, if is the Brownian loop measure (denoted by i/°°p in 
conditionally on the position of the SLEs, the probability that no loop in Lj intersects Ui<j7' and 7^ , 
7 = 2, . . . , 71 is: 

exp l^-A. ^({'5 ■ ^ ^ (U»<j7') ^ 0, <5 n 7^ 0}) j . 
It readily appears that the sum in the exponent can be written as: 

Y E in-lH{S:Sn-f^ ^0iS J e{ii,...,in}}) 

n>0 ii<---<in 

by an inclusion-exclusion argument. So it is equivalent to condition on no loop in Lj intersecting (at least) 
j distinct SLEs, j = 2, . . . ,n. 

Let us now formulate precisely and prove these results. Let k G (0, 8/3]. In the domain (H, xi, yi, . . . , x„, y„), 
consider n independent SLE^'s from Xi to yi, i = 1, ... ,n, and independent loop-soups L2, . . . , L„ with in- 
tensity Ak- The SLEs are defined by conformal equivalences {gl) and traces 7*. Assuming that the pairing 
{xi,yi), . . . , {xn,yn) is non-crossing, denote by iJ,{xi,...,yn) the law of the SLEs conditionally on the event 
that no two SLEs intersect and no loop in Lj intersects 7^ and Ui<j7* (if k < 8/3, one can drop the first 
condition). A hull A is a compact subset of EI such that H \ A is simply connected, A n R C A n H and 
xi,. . . ,yn are not in A] (pA is a conformal equivalence H \ A ^ A. Let L be yet another independent 
loop-soup with intensity A^; the union of A and loops that intersect it is denoted by A^; L^ is the collection 
of loops in L that do not intersect A. 

Proposition 1. Under the previous assumptions: 

(i) (Markov property) Under /ij-^.^ j,^), the law 0/ 17^(7^, . . . , 7j*j ^^-j, . . . 7") is that of (7^,..., 7") under 

M(3j(a;i)....gj(7j') gi{v„)) ("^P ^'^ time reparameterization) . 

(a) Let ip(xi, . . . ,yn) be the probability that no two SLEs intersect and no loop in Lj intersects 7-' and 
(U,;<j7*), j ~ 2, . . . ,n. Let gt = gl for some i G {1, . . . , n}. Then the process: 

i,{gt{x,), . . . ,g*(7D, • . •5t(yn))n U(^.)ff^fe) ( ^ ^ 

j^\ Katiy,) - 9t{xj) 

is a martingale under the (unconditional chordal) SLE measure for the i-th SLE , where = (6 — k)/2k. 
(Hi) (Restriction property) Under tJ,{xi,...,y„), the probability that A^ does not intersect (U7*) is given by: 

^(xr,...,y„) 1/ 1,^ ^"^^^ ^''^ U(x.)-0(y.) 

where (j) = (f>A- Moreover, the image under (j) of (71, . . . , 7„, L^, . . . L^, L"^) conditionally on "no loop in 
Lj connects 7^ to Ui<j7i, j = 2,...,n, and no SLE is connected to A by a loop in L" is distributed as 
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(7i, . . . , 7n, i2, ■ • ■ 7 where 7* are independent SLEs connecting 4>{xi) to 4>{xj), L2, ■ ■ ■ , L are independent 
loop-soups, conditionally on "no loop in Lj connects 7j to Ui<j7i, j = 2, . . . ,n". 

Proof. If 71 = 1, (i) is the usual Markov property for SLE, (iii) is a result from and (ii) is an empty 
statement (since ■(/'(a^ii yi) = !)■ 

By induction, assume that the results are proved for (n — 1) SLEs. Let us begin with (iii). Consider the 
event "no loop in Lj connects 7^ to Ui<:jji, j = 2, . . . ,n, and no SLE is connected to A by a loop in L". 
The probability of this event, conditionally on the n independent SLEs is: 

exp l^-x^ 1^ iy{{S ■.j^ns^0, {u^<jY) r\S^0}) + ^{{6 (u^j^) ns 0,An6 ^ 0}) 

For illustration, let us momentarily assume that n = 2. Then the sum in the exponential density reduces to: 

v{{6 : (5 n 7^ 7^ 0, (5 n 7^ 7^ 0}) + v{{6 : 5 n (7^ U 7^) 7^ 0, (5 n A 7^ 0}) 
which we can rewrite by inclusion-exclusion as: 

i^({(5 :5n7^7^0,(5nA7^ 0}) + v{{S : Sn-/^ ^ 0,SnA^ 0}) + iy{{S : 5n-f^^0,5n-f^^0,SnA^ 0}). 

The first term corresponds to the density of chordal SLE in (H\ A, xi, yi) w.r.t. chordal SLE in (H, xi, j/i); 
symmetrically, the second term is the density of chordal SLE in (IHI\A, X2,y2) w.r.t. chordal SLE in (H, X2, 2/2)- 
This follows from the restriction property for a single SLE, applied twice. The remaining term is the mass 
of loops in H \ A that do not intersect the two SLEs. So (iii) follows when n = 2. Let us get back to the 
general case: n > 2, and we assume that the assertions hold for (n — 1) SLEs. 

The sum in the exponential can be written as: 

^K{^:7^n<5 7^0,(U,<,Y)n<5 76 0}) 

j<n 

+iyi{S : 7" n (5 7^ 0, (Uj<„7') nS ^ 0,AnS ^ 0}) 
+2iy{{6 : 7" n (5 7^ 0, (U,<„70 D S ^ 0, An S ^ 0}) 
+i/({(5 : 7" n 5 7^ 0, (Uj<„7') n5 = 0,AnS ^ 0}) 
+iy{{S : 7" n 5 = 0, nS^0,An5^0}) 

and after rearranging: 

^ ^{{6 :j^nS^0, (U,<,y) nS^0}) + u{{6 : (U<„7^) nS^0,AnS^0}) 

j<n 

+iy{{6 : 7" n 5 7^ 0, ^ n (5 7^ 0}) 

+i/({(5 : 7" n (5 7^ 0, (Ui<„7*) nd^0,AnS^0}) 

The first line is treated by the induction hypothesis: it corresponds to the (unnormalized) density of 

H(m\A,xi w.r.t. /i(H,xi j/„-i) The second line corresponds to the density of 7" as an SLE in 

(H \ A,Xn,yn) W.r.t. 7" as an SLE in the full domain (H, a::„,j/„) (restriction property for a single SLE). 
The last line is the mass of loops in H \ A that intersect both 7„ and Ui<„7' (from the restriction property 
of the loop-soup). So the exponential above is indeed the density of fJ'{M\A,xi,...,y„) w.r.t. ^(H,xi,...,y„)- The 
statement on the joint law with loop-soups follows from the restriction property for the loop-soups L2, ■ ■ ■ ,L. 
(Note that we have exchanged loops between loop-soups, but these were only loops intersecting A). 
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The formula in (iii) is obtained by keeping track of masses of (unnormalized) measures in the induction 
scheme above; the new covariance factor in the product comes from the restriction property for a single SLE 
(say 7"). 

If A, B arc two hulls and A.B is the hull satisfying (pA.B = 0_b ° '/'a, then it is clear from the formula that 
the probability that Ui4>A{Y) does not touch B^, conditionally on "no loop in Lj connects to Ui<j7i, 
j = 2, . . . ,n, and no SLE is connected to A by a loop in L" multiphed by the probabihty that Lii(f>A{Y) does 
not touch is the probability that UiY does not touch (A.B)^ . This gives the (usual) restriction property 
for (Ui7*)^, which is less precise than the statement (iii). (See also |H] for a discussion of these restriction 
measures). 

Let us now prove (i),(ii). As explained above, the conditioning is in fact symmetric in the n SLEs. So, 
without loss of generality, we can assume that i ^ n. Reasoning as above, using the Markov property for 7" 
and the restriction property (iii) for (7^, . . . ,7""^) (where A is replaced with 7[J ^j), one gets easily (i) and 
(ii). ' □ 

When K = 2, this situation connects with Fomin's determinantal formulae f |12[ll5) '). Indeed, adding the 
loops of a loop-soup with intensity A2 to an SLE2, one gets a restriction measure with exponent 1, that has 
the same outer boundary as the Brownian Excursion. (In fact, the union of SLE2 with the unfilled loops 
that intersect it has the same distribution as the range of the Brownian excursion, see [H2IE1)- For instance, 
consider the situation with two SLE2's. Then the probability that no loop intersects the two SLEs is the 
same as the probability that the first SLE does not intersect the second one with loops attached, which is 
the same as the probability that an SLE2 does not intersect a Brownian excursion (with respective endpoints 
xi, X2, a;3, 3:4). This gives an interpretation of the symmetry of Fomin's formulae in the scaling limit (this 
symmetry follows from Wilson's algorithm in the discrete case). 

Under mild regularity assumptions, Ito's formula and (ii) imply that 1/', which is conformally invariant 
by construction, is annihilated by the n operators: 

-dl +y ^ + > + — + > 

Under the assumption of reversibility (see |23| ). ip is also annihilated by the n operators obtained by swapping 
the X and y variables. Note that the pairing is materialized in the constant terms of these operators. In fact, 
one can symmetrize the equations by a conjugation. Indeed, denoting I'^+j = j/j, the function 



V'(a;i,...,a;2„) J|( 



1-6/k 

\-7l-t-J -J/ 
J 



is annihilated by the operators 



7:dkk+y + > 7 r^, fc = l,...,2n 

2 ^-^ Xl ~ Xk K ^-^ [Xl ~ Xk) 

l^k l^k ^ ' 



Ek^k 

J2k ^kdk - n{l - 6/k) 

^ Y^k^k^k- {'^-&/n){xiA 'tX2n) 

the three first-order operators representing conformal covariance. This is the system (12. If) . So under regu- 
larity and reversibility assumptions, each non-crossing pairing of the {2n) points (xi, . . . , X2n) gives rise to a 
solution of this system. 

In the case where 8/3 < k < 4, one can expect that one can still define a natural law on non-intersecting 
SLEs, whose density w.r.t. independent chordal SLEs is given by exp(— AkI'(. . . )); this density is now 
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unbounded: 

d^J.{xl....,y„){'yl, • ■ • ,7n) =il^{xi, . . . ,y„)"^l{^,n7^=0,i<j} exp |^-A^^i^({(5 : 6 (Ui<j7') ^ 0,Sr\-y^ ^ 0})^ 

Xn,yn) iln) 

where dfi(^^. y.^ is the measure induced on paths (say as elements of the Hausdorff space) by standard chordal 
SLEk- So ip can be interpreted here as a (hopefully finite) partition function. 



4 Some particular cases 

Let K > and n > 1; xi, . . . , X2n are boundary points of H. Consider the operators (|2.1|l : 



2d, 



K — 6 



fc = 1, . . . ,2ri 



J2k ^kdk - n{l - 6/k) 
I Efc45fc-(l-6/t)(a;i 



If (/3 is annihilated by these operators and l is an involution of {1, . . . ,2n} with no fixed point (so that l 
determines a pairing), denote: 

Then t/i is conformally invariant (from the last three equations of (|2.1|l l: moreover, consider an SLE„ from Xj 
to 2^t(j), with associated conformal equivalences {gt) (with hydrodynamic normalization) and driving process 
Wt. Then: 



ip{gt{xi),...,Wt,...,gt{x2n)) n [9't{xk)g't{x,(k)) 

{kAk)]^{jAm V 



Xk - X,(k) 



gt{xk) - gt{x,(k)) 



is a local martingale (from the j-tli equation of (|2.1|) '). 

Wc study this system in a few cases, before giving a set of formal solutions. 



4.1 Case n = 2 

If tp and Ip are as above (say for the pairing {(xi,X2), (2:3, 0:4)}), then -0 is a conformally invariant function 
of four boundary points, so ip^xi, . . . , 2:4) = /(r), where r is the cross-ratio: 

^ _ {Xs - X2){X4 - Xi) 
{X3 - Xi){x4 - X2) ' 

Now it is easy to check that is a solution of (|2.1|l iff f{r)r^'^/'^ is a solution of the hypergeometric equation 
with parameters (a, 6, c) = (4/k, 1 — 4/k, 8/k) (see |1L)|). 

Assume that xi, . . . ,Xi are in cyclical order. It is easy to see that the boundary conditions ^^{xi, 2:3, xa, X4) = 
(r = 0), 'ip{xi^xi,xz,Xi) = 1 (r = 1) determine the following solution: 



r(4/At)r(i2A-i) (A 4 

1p[Xi,X2,X^,XA) = , , -T-'' ' 2i^l -,1 , 

I (8/k)I (8/k — 1) \n K 



4 8 
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for K < 8. For k = 6, one recovers Cardy's formula. For k S (0, 8/3], this can be interpreted as the probabihty 
that two slew's connecting xi,X2 and X3,X4 resp. do not intersect and are not connected by a loop in an 
independent loop-soup with intensity A^- 

A solution of the hypergeometric equation with parameters (a, b, c) = (4/k, 1 — 4/k, 8/k) can be written 

as: 

' U-4/«(l_^)12/«-2(l_^2.)-4/K^^^^l-8/K r „-4/«(„_.^)12/K-2(^_^)-4/«^^ 

c Jc 
where C is a formal linear combination of "cycles" , i.e. paths of integrations over which the integrand is 
single-valued and are either closed or start and end at singular points of the integrand, i.e. (0, l,x,oo). In 
terms of the original variables, one gets the following expression (inverting the roles of x^ and X4): 

^(Xi, X2, X3, X4) = n (^^- - ^^)'^" 11(^4 - X,y~'^^ j \{{V - X,)-'/^V - X^f^'^'^dv (4.2) 

Note that although (|2.1|l is symmetric in the (2n) variables xi, . . . ,X2m the integrand is symmetric only in 
the first three variables. In fact, one can interchange the role of X4 and, say, X3 by a homographic change 
of variables. It is possible, but apparently not very practical, to rewrite this integrand in a completely 
symmetric fashion. 

4.2 Case k = 2 

If one specializes the above expression to k = 2, one gets a rational function, since 2^1(2, —1, 4, r) = 1 — r/2; 
note that 2r(l — r/2) = 1 — (1 — r)^. The corresponding solution ip is: 

(p(Xi,X2,X3,Xi) 



{Xi - X2Y{xz - 2-4)2 (xi - X-iY(x2 - Xi)^ 

(a;i-a;2)~^ (xi - x^y'- 
ix4-X2y^ (xi-xsy- 



det 



For a general value of n, it is possible (if a bit tedious) to check that: 

(p{xi, . . .,X2n) = det {{Xi - Xj+ny^)^^^ ,^^^ 

gives a solution of H2.1|l . corresponding to the pairing {(cci, a;„-)-i), . . . , (x„, X2n)}- Other pairings gives 
different solutions; there are non-trivial linear relations between those solutions. These solutions correspond 
to the scaling limit of Fomin's formulae ( |12II15| ). More precisely, multiplying the determinant by the product 
of its diagonal terms, one gets an alternating sum of probabilities of a certain non-intersection event. It is 
interesting to observe that if {xi , Xj ) are paired, the limit of this probability as Xi — *■ Xj is given by a lower 
dimensional determinant, as it should. Likewise, the decorrelation of SLEs living on different scales is obvious 
from these formulae. For the particular pairing {xi, X2n), ■ ■ ■ , (xn, Xn+i), if the points are in cyclic order, the 
probability of the corresponding non-intersection event is given by a single determinant. 

Each permutation of xi, . . . , X2n in this formula give a solution of the system; there are (many) non-trivial 
linear relations between those solutions. In general, other (geometric) pairings do not correspond to a single 
determinant, but to a linear combination of determinants of this type. One has to determine the solution 
satisiying appropriate boundary conditions. It is clear that each determinant (divided by its 

diagonal terms) goes either to or to a lesser order determinant not involving Xi,Xj. 
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4.3 Case k ^ oo 



Consider the following system: 




k = 



n 



1, 



,2n 



+ X2n) 



which can be seen as the (somewhat degenerate) limit of (|2.1|) as k goes to infinity. Then it is easy to see that 
solutions (functions annihilated by those operators) are polynomials in xi, . . . X2n homogeneous of degree n, 
of partial degree at most 1 in each variable. The vector space of these polynomials has dimension C(2n, n); 
the dimension of the subspace of such polynomials that are also translation invariant is 



Considering the derivative of these polynomials with respect to a variable, and by induction on the number 
of variables, one gets that this linear space is spanned by: 



and polynomials obtained by action of the symmetric group on this one. Note that these spanning polyno- 
mials also satisfy the last relation in the system. So the solution space has dimension exactly C„. 



Here we consider the system satisfied by the crossing probabilities for critical percolation in an hexagon with 
alternate boundary conditions (or a simply connected domain with six points marked on the boundary). In 
this situation, there are C3 = 5 elementary crossing events, with probabilities adding up to 1 (see Figurc[5J|. 
A solution (p of this system is conformally invariant (and not only covariant); so one can send xi,X2,xq to 
00,0, 1 by an homography, and tp becomes a function of three real variables 2/i,?/2,?/3- This solution must 
be annihilated by six second-order operators; it turns out that there is (exactly) one linear relation between 
these operators. Singularities for these equations occur when two of the y variables are equal, or one is equal 
to 0, 1, 00. Consider an open set U of regular points. Let D{U) be the subalgebra of differential operators in 
C(2/i, ?/2, y3)[di,d2, 83] with coefficients regular in U, and / be the left ideal of differential operators generated 
by these equations. 

It is possible to write the system as an integrable, rank 5 Pfaffian system (see e.g. this is also true 
for any k > 0). Also, note that the constant function is a solution; we will use this simple fact to explicitate 
solutions of the system. Consider the image of the the solution space in U (one may also think in terms of 
local systems, sheaves of algebras, ...). Its image under the operator di has dimension 4. Let J be the left 
ideal: 



Then J is a holonomic ideal of rank 4. There is a well-known example of a holonomic system of rank 4, with 
3 variables, and the same singular locus as J, namely the system of Lauricella's Fjj function, which is the 
natural generalization of Appell's Fi (see e.g. |1(JI IT]'). 

Recall that the system Fd with parameters a,7,/3i, . . ./3m in the variables wi, . . . ,Mm, is given by the 
linear operators: 



C(2n,n) - C(2n,n- 1) = C{2n,n)/ {n + I) = Cn- 



(xi - Xn+l 



)■■■ {Xn - X2n) 



4.4 Case k = 6, 2n = 6 



J = {Le D{U) : Ldi e /}. 



m 



Uj{l - Ui)dij + (7 - (a + /3i + l)ii.i)9j - ^ I3iu-jdj - a(3i 
for i = 1, . . . , m. A solution then also satisfies the Euler-Darboux equations (if 7 — a 



1^0): 
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for I < i < j < m. It is known that this system is of rank (m + 1) and that an integral representation of 
sohitions is given by: 



By ehmination in D{U), it is possible to exhibit elements in J that generate a rank 4 ideal, the solution 
space of which is a conjugate of that of a particular system. More precisely, with the help of a formal 
computation software, one can prove that if (ys is a solution of the original system, then: 

(yi(2/2-2/i)(y3-2/i)(i-yi))^/% . ^ 



((1 - y2)(i - y3){y3 - 2/2)^22/3)^/^ 

is a solution of the Fo system with parameters: 

a = 1/3, 7 = 2/3, A = -4/3, P2 = 2/3, P3 = 2/3. 

Similarly, in order to get a solution symmetric in the y variables, one can do the same computation for the 
operator {yidi + ?/2C?2 + 2/3^3), which by conformal invariance corresponds to a perturbation of xq. Then it 
turns out that: 

((l-yi)(l-2;2)(l-2/3))'/' ^ «^ ow ^ 

(yidi + y2d2 + y3d3)nyi,y2,y3) 



{yiy2y3{y2 - yi){y3 - yi){y3 - y2)Y''^ 

is a solution of the Fo system with parameters: 

a = 1/3, 7 = 8/3, /3i - 2/3, (32 - 2/3, P3 = 2/3. 
Solutions of this system are given by 

u-2/3(l - uf/^ ((1 - uy{){l - uy2){l - ^/y3))"'/' 



where as before C is an adequate "cycle" . To get to we write: 

-^ip{syi,sy2,sy3) = {{yidi + y2d2 + y3d3)ip)isyi, sy2, sy3) 
which after trivial manipulations leads to the following expression for ip{xi, . . . ,xq): 

n (^^- - ^•')'^' 11(^6 - x,f j \{ {u^- x.r^/s -Q _ ^^)0(^^^ _ ^^f/^du^du2. (4.3) 

4<i<6 i<6 '-^ie{l,2} ie{l,2} 

We end this section by quoting a particular result for configurations with symmetries. Consider the 
domains {lJ,l,u,j,ju,j^,j^u) where U is the unit disk, j = e^"/-^, and u is on the arc (l,j). Restricting a 
solution Lp to such domains, one gets a function g of a single variable u. Now, computing in D(U) as above, 
it is possible to derive a third-order ODE satisfied by g (of course this ODE has no constant term). The 
fact that the rank goes down is linked to the particular symmetries of the configurations we study, and is 
easily interpreted in terms of percolation obscrvablcs (different elementary events have the same probability 
because of the threefold rotational symmetry of the configuration). Making the cubic change of variables 
V ~ (that sends the singularities to 1), then the quadratic change w = — (w — 1)^/411 (that sends 

the singularities 0,oo to 00), so that w = sin^(30/2) if u = e'^, it turns out that w f-* w^^^g'{w) satisfies 
the (classical) hypergeometric equation with parameters (a, 6; c) ~ (5/6, 5/6; 7/6). Hence g' belongs to the 
vector space spanned by: 



hi{w) = w' 
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Let gi{'w) ~ hi{s)ds. The function g belongs to the vector space spanned by (1, 171, (72)- Let us work out 
boundary conditions for crossing probabiHties. Consider: 

the probabihty that the three blue sides belong to the same blue cluster (the "Mercedes" configuration) . Let 

9+{w) = g{u) + g{j/u),g^{w) = g{j/u)~g{u) 

for u in the arc (l,e"^/'^) (note that when u goes from 1 to j, w goes from to 1 and then back to 0). 
For parity reasons, it is easy to see that one can write g+ = cigi + C3, g- = 02(72- We need to determine 
the three constants ci, 02,03. They are fixed by the boundary conditions: gj^{Q'^) = (7_(0+) = 1, and 
g+{w) — g-{w) = 0{y/w) = 0{u — 1) as w \ (the exponent here is three times the half-plane one-arm 
exponent). From the Euler integral: 

3^2(01,02, as; PiP2;x) = B(ai,pi-ai)-iB(a2,P2-a2)"' / s''^-Hl^s)P^-''^-H''--\l-t)P--''^-Hl-stx)-''^d.sdt 

J[0,l]2 

one gets: 

gi(0) = 5(5/6, 2/3)-i f\l - s)-^/^ f ^-^/^(l - - sty^^^dsdt 

Jq Jq 

= 3i^2(l, 5/6, 5/6; 3/2, 3/2; 1)5(1, 1/2) 

32(0) = 5(1/3, l/6)-i / (s(l-s))-i/2 / t-^/^i-t)-^/(i(l-ts)-^^'^dsdt 
Jo Jo 

= 3^2(1/2, 1/3, 1/3; 1, 1/2; 1)5(1/2, 1/2) = 2^^(1/3, 1/3; 1; 1)5(1/2, 1/2) = ^^^^^^^i'^^' ■ 
For the condition (ci/ii — C2^2)(w) = 0(w^^/^), we need the analytic continuation formulae: 

1/2 r(3/2)r(-i/6) p /5 5 7 \ , r(3/2)r(i/6) ^(225 



^ r(2/3)2 V6' 6' 6'"^; ' r(5/6)2 '^'V3'3'6 

1/2, , X r(i/2)r(-i/6) ^ /5 5 7 \ r(i/2)r(i/6) ,,6^/225 \ 
" ' '^^("^ = r(i/6)^ (e ' 6^ 6 ^ " j + r(i/3)^ ^ (3 ' 3 ^ 6^ V 



Hence we can conclude that: 

r(2/3)^ V3^,„,„,3 2r(5/6)^ / V3 \\,^,^,9 

^ r(l/3)r(l/2)^ - 2;^^^2/3) ,ci C2^;^ = j r(2/3) . 

So the probability of a configuration where two (given) blue sides and two yellow sides are connected for 
a regular hexagon {u = e^^^^, w = 1, see Figure ISJ is: 



I-C3 ^ ci.gi(0) ^ 2 / V3 
3 3 3 i 22/37r 



r(2/3)V2(l,5/6,5/6;3/2,3/2;l). 



5 Formal solutions 

In this section, we discuss a family of formal solutions for the system (|2.1f) . These are integrals taken on 
a "cycle" ; they define actual solutions under conditions of convergence and determination of the integrand. 
We shall discuss later the problem of identifying probabilistic solutions. 
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Figure 5: Percolation in a regular hexagon 



Consider the function: 

l<i<n l<i<n l<i-i<i2<n 

l<j<2n 

n i^h-^h)^ n (2;2n - 

l<Jl<j2<2n l<j<2n 

where x = {xi, . . . , a;2n), u = (wi, . . . We want to prove that a function of the parameters (xi, . . . X2„) 

of the form: 



/ (/)(x, u)du 
Jc 



defines a solution of H2.1|l for appropriate cycles C (that may depend on the x parameters). Note that this 
covers the cases (|4.2|) and (|4.3|) . The conformal covariance conditions are easy to check. Denote by the 
differential operator: 

'"2 ^xz-xfc K ;^(-^'-^^0'' 

Then the following lemma holds. 
Lemma 2. (i) If 1 < k < 2n, then: 

(ii) Moreover: 



dui \u.i - Xk 



\ ^ d I I Z, ft — O T-r U," — Xi -I r 

/:2n0 = - E 7^ — + — n — — - n 



X2n - Uj 



20 



Proof, (i) By symmetry, one can assume that fc = 1. Then: 



Cicj) K 



E 



E 



l<i<7i l<i<2n ^ 



1 - 



K / Xi — X2r. 



4 >^ 1 2 1 



6 



- 1-- 



1 



E 

l<j<2n 



X2n - Xi 



E 



H J {xi - X2n)'^ 

6 



k^3 



1-^ 

K 



- E — 

K ^ — ' Xi ~ 

l<k 



KKn l<j<2n 



which after simphfications leads to: 



2 1 



^ 1 8^/1 1 



l<i<n 



l<j:<ri 
l<j<2n 



1 



\<j<2n ^ J 



1 



l<i<n 



1 



1 



1 1 



16 



E 



Xl — X2n Xl — Ui X2n — Xi X2n — Ui 
1 



K Xl Uij Xl 

l<ii<i2<ri 



where we have used the identity: 



1 



(a-6)(a-c) {b-a){b-c) (c-a)(c-6) 
for (a, 6, c) = (xi,Xi,Xj). 1 < i < j < 2n. Also: 



1 a / 20 



9ui \Ui — Xl I Ui — Xl 



1^-2 



- 1 



1 



' ' ^ \<j<2n ' J 



1 8 1 
+ - > 

I<i2<n 
i2/i / 



E 

:2n 
\ 



Making use of the same identity as above, one can easily conclude. 

(ii) The situation here is slightly more intricate than in (i). Let P be the rational function: 



P=-\C 



2n(l 



E 

l<i<n 



d f 2(1) 



dui \Ui- X2n 
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Reasoning as above, one gets the following identity: 



V - 2 V ^ 



{u^-X2n)^ 4^ (U^ - X2n){Xj - X2n) 

l<j<2n 



y t + 4 y t 



On the other hand, if Q is the rational function: 



l<j<2n ^" ^ i<j<n,j^i ^ ■' 



2 

X2„ - Uj 



then Q can be written as: 



^ _ (8 - k)(4- k) ^ 1 -Q M» - J-j- (x2n~Uj'^'^ 

l<i<n l<j<2n l<j<7i,j^i 



E ^— + ^— + 2 V 



— Xi Uj — X2n — Ui — Ui 

So the statement reduces to the identity of rational functions P = Q. This is a bit tedious, and we include 
these computations for the sake of completeness. We can assume that k ^ 4,8. Note also that P and Q are 
symmetric in the u variables. Let us expand Q at ui = X2n (denote e ^ ui ~ X2n)'- 



K. P \ ^ ^ T"o„ — T" ■ ^ ^ 



l<i<n ^ ' , 



1 1 + ^ y + f y 7 ^ T + y 

\ l<j<ri \l<ii<i2<n ^' l<i< 

E — ^ — E — 

\ KKn l<j<2n 

3 1 / 3-1 2-6 
> + > 



l<j<2ri ■> \<i<n 

E3 — 2 ^ — A 12 
/'to.. T . If To.. T .- 1 ^- ^ 



l<J-l<j2<2n " ^Ji)(^2„ " X ^^^f^^^^^ {x2n " tin)(^2n " "i, ) 

v^9-2-4 1-1 -2 + 2-2 

+ / 7 T9+ / 7 T9+ / 7 77 7+o(e). 

,7~^ (a;2n-Mjj i/"^.-, \X2n-XjY ^ \X2n - X j)(X2n - Ui) 

l<j<2n 

From here, and using symmetry, it appears readily that: 

P -Q = o{{ui - X2„)) 



22 



for 1 < i < n. One can also check that Q is regular along u,; = Uj, for 1 < i < j < n (by symmetry it is 
enough to check that the term in {ui — Uj)~'^ vanishes). Considering {P — Q) a,s a. rational function of the 
u variables, it appears that it has no pole (even at infinity), so it is constant (i.e. a function only of the x 
variables). Setting ui = X2n, it appears that {P — Q) is identically zero, which concludes the proof. □ 

One can rephrase the lemma as follows: there are differential (n — 2)-forms wi, . . . ,u)2n with rational 
coefficients such that for 1 < fc < 2n: 

Ck{(t)dui A • • • A rfM„-i) = d[cj)ujk)- 
So for a "cycle" C, we get a (real) solution of (|2.1() as soon as the following formal computation makes sense: 

Ck I (pdui A ■ ■ ■ A dun-i ^ I Ck{<l>dui A ■ ■ ■ A dun-i) = / d{(j>u!k) = / ^Wfc = 0. 

Jc Jc Jc JdC 

The cycle C must be such that cj) has a single-valued determination on C, one can differentiate w.r.t. the 
X parameters in the integral, and (f)UJk vanishes (to a sufficient order) on dC (which is true in particular if 
dC = 0). In fact C is a cycle for the twisted de Rham homology associated with (see e.g. section 
5.4). In the next section we give examples of such cycles, in relation with the probabilistic set-up. 



6 Explicit solutions 

In this section we propose choices of cycles of integration such that one gets well-defined solutions, and that 
these solutions can be interpreted in the probabilistic situations discussed earlier. We distinguish the cases 
K € (0, 8/3] and n — 6. Recall that (with an emphasis on the number of parameters): 

n— — TT ——2 TT — 

{Ui~Xj) " _[_[ {Ui-X2n)'' [[ {Ui^-U^^)" 

l<i<n l<i<n l<ii<i2<n 

l<i<2n 

n fe-a^ii)^ n {^^n-xjf~^ 

l<il<j2<2n l<j<2n 

where x = (a;i, . . . , X2n), u = (ui, . . . u„_i). 

Assume that xi < ■ ■ ■ < X2n are boundary points of H, and the involution t with no fixed points define a 
non-crossing pairing of the (2n) points xi . . . X2n- Recall that the fundamental group ni(C \ {xk, x^^k)}, zq) 
is the free group generated by two elements ak, o'tik) corresponding to loops around Xk, a^t(fc) respectively. A 
double contour loop around Xk,x^(^k) corresponds to the commutator 

(^k<y,{k)<yk^cr;-J^f 

so that its image in ni(C \ {.Xfc}, zq) and ni(C \ {^^(fc)}, zq) is the identity. We assume that the other x 
points are outside this loop. Then is single- valued on such a double contour loop. Now let us enumerate 
the (n — 1) pairs that do not contain the last point X2n' 

{{xk,x,{k)},k = 1, . . . ,2?!} = {{ak,bk},k= 1, ... (n - 1)} U {{.T2„, Xt(2„)}}. 

Denote also a„ = Xi,(^2n)j bn = X2n, and assume that Cfe < 6^, fc = 1, . . . n. For fc = 1, . . . , n — 1, let Cfc be a 
double contour loop around ak,bk, so that these loops do not intersect pairwise (see Figure E)). Define C to 
be the Cartesian product of these loops: 

C = Cl X • • • X Cn-l- 
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Figure 6: Two non- intersecting double contour loops 



Of course. C is a function of xi, . . . , X2n- Then 0„ (as a function of ui, . . . , Un-i) is single- valued on C, and 
dC = 0. Moreover, the integral of (/)„ on C does not change if the loops are deformed. 

If i is a double contour loop around 0, 1, and p, q are numbers with real part larger than (—1), then: 

T{p + q + 2) 

as is easily seen when L gets close to the unit segment. By analytic continuation in p and q, the second 
expression stays valid for general values of p and q. In the situation where p ~ q, one can use figure eight 
loops. 

Assume that the two consecutive points Xk,Xk+i, fc + 1 < 2n, are paired and correspond to the loop Cj. 
Then it is easy to see that: 

lim (xk+i - Xk)^^"^'^ / 4>„{:x.,u)dui . . .dun-i ^ c 4'n-i{x,u)dui . . . duj . . . dun-i 

where x = (xi, . . . ,Xk,Xk+i, . . . , X2n), u = (ui, and C = Ci x • • • Q x • • • x C„_i. The 
constant c is given by (up to a complex number of modulus one depending on the choice of a determination 
for (pn on C): 



c = Ck = 4 sin — TT 



4 \2r(l-l)' 47r2 



r(2-f) r(2-f)r(i)^ 

This constant is non zero, since we assume that (8/k) ^ N. If n = 2, then the limit is c{xi — xzY~^/'^. 
In this set-up, we can now formulate: 

Proposition 3. Let k € (0, 8/3), 8/k ^ N. Let 71, . . . , 7„ he (the traces of) n independent SLE^ 's from Uk 
to bk, k = 1, . . . n, and L2, ■ ■ ■ Ln be independent loop-soups with intensity A^. Then: 



n (U,<,7. = 0), J = 2, . . . ,n) = (cji- \[{bk a,f/--' 



fe=i 



c 



i„(x,u)du 
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Proof. Denote by "tp^xi , . . . X2n) the right-hand side of the equation in the proposition. Then we want to prove 
ip = ip (with the notations of Proposition^. By construction, ip has the property (ii) of Proposition^ Also, it 
is easy to sec that -0 is conformally invariant. So let us momentarily fix the values xi = 0, X2n-i = 1, X2n = oo. 
Then ijj extends continuously to the compactification of 

{{X2, . . . , X2n-2) e K^""^ S.t. = Xi < • • • < X2n-1 = 1} 

so it is bounded. There is a fc € {1, . . . , n— 1} such that and bk are two consecutive points on the real line. 
Consider the martingale associated with ip and the evolution of the fc-th SLE. We use the hydrodynamic 
normalization at infinity; (g^) arc the conformal equivalences associated with the SLE, which is defined up 
to time T — T{bk) (which is finite due to this arbitrary choice of normalization). Then g^^ defines a conformal 
equivalence between the unbounded connected component of H \ 7^. and H. As t r, Wt and gt{bk) go to 
the same finite limit, and gt{aj), gt{bj) have finite (distinct) limits. So the martingale goes to: 



^n-ligrmUig'rMg'ribj) 

with the obvious notations. The product can be written as: 



bj - a J 



9r{bj) ~ gr{aj) 




where L is an independent loop-soup. Now, as we remarked earlier, ^/'(x) does not depend on the ordering 
of the SLEs (this follows from the Poissonian nature of the loop soup). So we can assume that the [n — 1) 
SLEs loop-soups are reordered so that the last one corresponds to k (that is, fc is re- indexed as n). As in 
Proposition ^ the restriction property for the loop-soup and for (n ~ 1) SLEs now implies that: 

V'(x) = E P ((7,)^^ n (U,<,7.) = 0, J = 2 . . . (n - 2)\g^{a,),gr{bj),j < [] \^,)^^n7„ = 
\ i<" ' 

where 7^ is an SLE from grio-j) to g^ibj), and ij, L'^ are independent loop-soups. Comparing with the 
limiting value of the martingale as i /" r, one concludes by induction on n, using the optional stopping 
theorem. □ 



Note that it appears to be difficult to determine analytically the boundary behaviour of ip along all 
boundary components. The point here is that the dynamics of the SLEs lead to an exit on a fixed boundary 
component; and commutation allows us to consider the SLEs in an appropriate order. 

We gave probabilistic interpretations in the case n < 8/3, k = 6, k = 8. We have identified corresponding 
integration cycles when k < 8/3, 8/k 7^ N and k = 8. Let us discuss the remaining cases k = 8/m, m > 3, 
and K = 6. 

The problem is to identify a cycle of integration corresponding to a geometric configuration of SLEs. 
In particular, one can consider the case of n nested paths connecting Xi to X2n+i~i, i = 1, . . . ,n, where 
.Ti, . . . ,X2n are in cyclical order on the boundary. This is the natural situation for Fomin's formulae (but 
not for the UST e.g.). Then the following cycle can be used (say when k < 4): 

C = Ci X • • • X Cn-l 

where Ci is a loop starting and ending at X2m circling counterclockwise around X2n-i^ and leaving the other 
marked points on its right-hand side; moreover, the C^'s are chosen so that they do not intersect except at 
X2m consequently the integrand has a single- valued determination on C. 
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In the case where k = 2, induction on n and a residue computation give the determinant obtained as the 
scahng hmit of Fomin's formula. One could also have chosen C = C[ x ■ ■ ■ x C'„_i, with Ci a loop starting at 
xi and circling around xi^i (and the special point in the integrand is now xi, not X2n)', the solution is then 
proportional to the one corresponding to C. Here, using induction on n (and analytic continuation in k), we 
can show that these two solutions are indeed proportional and have the correct asymptotic behaviour when 
two consecutive marked points collapse. Finding a cycle with geometrically prescribed boundary conditions 
(in the Weyl chamber {xi < ■ ■ ■ < X2n}) seems to be technically difficult in general. 

Also, we remark here that the integral representation becomes somewhat degenerate in the particularly 
interesting case k = 6. As in the case of standard hypergeometric functions, analytic continuation in the 
parameters (here k) can be put to good use. Since the functions ip can in this instance be interpreted as 
(scaling limits of) crossing probabilities, it would be nice to find real integration cycles (formal sums of 
Cartesian products of segments [xi,Xi+i]) corresponding to these events (i.e. satisfying the right boundary 
conditions; in general, when two consecutive points collapse, the boundary condition is either or a crossing 
probability for a configuration with 2(n — 1) marked points). For small values of n (e.g. 2n = 6), one can 
work out such cycles, though the general pattern is not very clear. Again for small n, these cycles give 
enough solutions, so the results of this paper can be seen as giving an algorithmic solution to the problem 
of computing crossing probabilities for the scaling limit of percolation with alternate boundary conditions. 

As for general properties of the holonomic system we studied, the main question is to prove that the 
rank is indeed C„. From the case k — > oo, we see that the rank should be at most Cn- Under smoothness 
assumptions, the restriction construction give C„ distinct solutions when < k < 8/3; this should be enough 
to prove that the rank is C„ in general (since all these systems then have a rank Cn Pfaffian form that should 
have an analytic continuation in k). 

Acknowledgments. I wish to thank David Wilson and Greg Lawler for fruitful conversations, as well as 
anonymous referees for helpful comments. 
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